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I. INTRODUCTION 

In contrast to point contacts between normal metals, 
the superconducting point contacts have very non-trivial 
transport properties even under the simplest bias condi- 
tions of fixed dc bias voltage V across the contact. The 
origin of this complexity is the oscillating Josephson cur- 
rent induced in the point contact by non-vanishing bias 
voltage which makes electron motion through the con- 
tact inelastic. From the microscopic point of view, this 
type of inelastic electron dynamics can be understood 
in terms of cycles of Andreev reflections at the two in- 
terfaces between the contact and superconducting elec- 
trodes Q|. Electron emerging from, say, the left elec- 
trode with energy e relative to the Fermi level of this 
electrode is accelerated by the applied voltage V and 
reaches the right electrode having the energy e + eV. 
It is then Andreev-reflected back as a hole which is again 
accelerated when crossing the contact backward. After 
being Andreev-reflected from the left electrode this hole 
produces an electron with the energy e + 2eV (Fig. la). 
Because of this process of multiple Andreev reflections 
(MAR) electron with energy e incident on the contact 
can absorb or emit some number n of quanta 2eV of the 
Josephson oscillations and emerge from the contact area 
with the energy e + 2neV. 

Phcnomcnologically, the cycles of MAR manifest them- 
selves in the so-called "subharmonic gap structure" 
(SGS) in the current-voltage characteristics of the su- 
perconducting contacts: current singularities at voltages 
V n = 2A/en, n — 2,3, where A is the superconducting 
energy gap of the contact electrodes. Since the cycle of 
two Andreev reflections requires that both electron and 
hole traverse the contact, the probability of this process 
is proportional to D 2 in contacts with small transparency 
D. This means that injunctions with the low-transparent 
barriers, the MAR current is much smaller than the 
regular quasiparticle current or current of Cooper pairs 
which both scale as D with junction transparency. There 
is, however, a very large variety of high-transparency 
Josephson junctions, where multiple Andreev reflections 
and associated subharmonic gap structure are important. 
Possible implementations of the high-transparency junc- 
tions include tunnel junctions with high critical current 
density @], semiconductor/superconductor heterostruc- 
tures (^J4|], controllable break junctions disored- 
ered superconductor/semiconductor junctions short 
metallic SNS junctions IlC(]. 
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FIG. 1. (a) Schematic energy diagram of the multiple An- 
dreev reflections process in a ballistic constriction between 
two superconductors with energy gap A. Both electrons and 
holes are accelerated by the applied bias voltage V and gain 
energy eV crossing the constriction, (b) Geometry of the 
constriction. The shaded area corresponds to the scattering 
region. The pair potential A(x) can be neglected in regions I 
and II. 

The subharmonic gap structure was extensively stud- 
ied in experiments done in the 60's Jll]-[l5| , and was inter- 
preted in terms of the multi-particle tunneling Q . How- 
ever, the theory of multi-particle tunneling is valid 
only in the limit of low junction transparency, D — > 0, 
and can not account quantitatively for all features of the 
SGS. A more detailed understanding of the SGS was de- 
veloped later |I|,|l9|-^l| based on the idea of MAR. Never- 
theless, quantitative dependence of the SGS on the con- 
tact transparency D was still not fully understood at that 
time. During the last few years, considerable progress has 
been made in quantitative description of MAR in short 
constrictions (shorter that the coherence length £ of the 
superconducting electrodes of the constriction) with arbi- 
trary transparency |2^-|26|] . The aim of this Chapter is to 
review these recent developments. Section 2 presents the 
scattering approach to MAR in the multi-mode constric- 
tions, in particular, in short disordered SNS junctions. 
In Section 3 we discuss the Green's function approach to 
description of MAR which is needed to account for the 
effects of non-trivial microscopic structure of the con- 
striction electrodes. In Section 4 the noise properties of 
the ballistic constrictions are calculated starting from the 
results obtained in Section 3. The Conclusion summa- 
rizes the results and unsolved problems of the theory of 
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ac Josephson effect in high-transparency Josephson junc- 
tions. 



II. MULTI-MODE CONSTRICTIONS 

The main simplification brought about by the con- 
dition d <C £, where d is the characteristic dimensions 
of the constriction, is the possibility of neglecting, first, 
all properties of the constriction besides its normal-state 
scattering matrix S, and second, all deviations from 
equilibrium in the electrodes of the constriction |^,|l9| . 
These two factors allow to avoid the self-consistent de- 
termination of both the pair potential and electrostatic 
potential. Electron transport through the constriction is 
determined then by the interplay of scattering in the con- 
striction (described by S) , acceleration of quasiparticles 
in the constriction by the bias voltage V, and Andreev 
reflection with the amplitude a(s) at the two interfaces 
with the superconducting electrodes in equilibrium, 



modes, all amplitudes of the electron and hole wavefunc- 
tions have transverse mode index m not shown in eqs. 
(g) and (§), e.g., A n = {A n , m }, m = 1, N. The source 
term J describes an electron generated in the jth trans- 
verse mode by a quasiparticle incident on the constriction 
from the left superconductor: J(e) = (1— \oq | 2 ) 1//2 <5 mj -. 

The current in the constriction can be calculated in 
terms of the electron and hole wavefunctions in the region 
I or II. The contribution i(t) to the current from the 
wavefunction (@) is: 



*(*) = — ImTr(VezV^ - VftWl) , 



(4) 



where Tr is taken over the transverse modes. Equa- 
tions (||) and (|4|) imply that the total current I(t) in 
the constriction oscillates with the Josephson frequency 
ujj = 2eV/h and can be expanded in the Fourier compo- 
nents: 



ikuj i 



a(e) = 



e-sign(e)(e 2 - A 2 ) 1 / 2 . 
e-i{A 2 



2\l/2 



£ |> A. 
£ |< A. 



(1) 



Here e is the quasiparticle energy relative to the Fermi 
level of the electrode, and A is the equilibrium energy 
gap of the electrodes. 

To find the transport properties of the constriction 
quantitatively, we generalize 1 23 2^] to non- vanishing bias 
voltages V the scattering approach for Bogolyubov-de 
Gennes equations |2|-|3l|] that describes dc Josephson ef- 
fect at V — 0. In accordance with the qualitative picture 
of MAR discussed in the Introduction, electron with en- 
ergy £ incident on the constriction generates electron and 
hole states at energies £ + 2neV with arbitrary n. Thus, 
the electron and hole wavefunctions in regions I and II of 
the constriction (Fig. lb) can be written as follows: 



Substituting eq. (Q) into (Q) and summing the contribu- 
tions from quasiparticles incident both from the left and 
right superconductors at different energies £ we obtain 
the Fourier components Ik of the current: 



h = 



e 

itTi 



-eV 



detanh{— }Tr[(Jjt4 - 



'2k 



(1 + a 2n a* 2(n+k) )(A n Ai +k - B n B] l+k ))\ 



(5) 
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ik.t: 



^ ~Q g — ikx^ — i(s J r2neV)t/h 
l e -i(e+2neV)t/h 



The amplitudes A, B, C, D of electron and hole wave- 
functions (|^) and (^) are related by the matrix S of scat- 
tering in the constriction. Taking into account that the 
, scattering matrix for the holes is the time-reversal con- 
jugate of electron scattering matrix 5* we can write: 



(2) 



(II) 



c ikx 



-ikx^ e -i(e+(2n+l)eV)t/h 
-ikx] -i(c+(2n+l)eV)t /h 



(3) 



B n 

C n 

A n 

Dn-l 



= S 



(l2nA n + J S n o 
0>2n+lDn 



S* 



a 2n B n 



(6) 



(7) 



where a n = a(e + neV). In these equations we took 
into account that the amplitudes of electron and hole 
wavefunctions are related by the Andreev reflection, and 
also neglected variations of the quasiparticle momentum 
k with energy assuming that the Fermi energy in the elec- 
trodes is much larger than A. The quasiparticle energies 
in regions I and II are measured relative to the Fermi level 
in the left and right electrode, respectively. Since the con- 
striction is assumed to support TV propagating transverse 



The scattering matrix 5* is a unitary and symmetric ma- 
trix 2N x 2N and can be written in terms of reflection 
and transmission N x N matrices r, t: 



S e i = 



r t 
t' r' 



(8) 



where t 1 — t T . 
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Eliminating A n between eq. (|6| 
we find the relation between the amplitudes B n and D n 



rr' =1. 
and inverse of eq. (|t]) 
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Combining this relation with the expression for D n in 
terms of B n that follows from the inverse of eq. (||) and 
eq. (0) we arrive at the following recurrence relation for 

B.„ .: 



tit 



B n +1 — (" 



4 2n+l 



''2d 



*2n+l 



ft: 



2n+l 



£? n _i ) - [1 - a\ n \B n = -rJS n0 , (9) 



a2nfl2n-l 

+ 1 _ n 2 
1 a 2n-l 

In a similar way we obtain the recurrence relation for the 
amplitudes A n : 



,»,T / a2n+l02n , 
' ' I 2 — 

1 — a 2n+l 

a2n-l02n-2 



2 ft 2 
,„ + , °2n-l )An , 



t °2rt-l 



t a 2n+l 1 a 2n-l 



(ai(5„i - a <5„ ) J 



1 - a 2 
-aoJS n o ■ 



(10) 



Since the hermitian matrix iit can always be diago- 
nalized by an appropriate unitary transformation U, the 
recurrence relation (^|) implies that the structure of the 
amplitudes B n as vectors in the transverse-mode space 
is: 



B„ 



tff n (D)UrJ, 



(11) 



where D = U ttfU^ is the diagonal matrix of transmission 
probabilities D m , m = 1, N. The functions f n {D) are 
determined by the solution of the recurrence relation (||) 
with the diagonalized transmission matrix itt. 

Equation (|ll|) shows that the contribution of the am- 
plitudes B n to the currents (||) can be written as 

Tr[B n Bl] = (1- \a \ 2 )Tr[f n (D)f*,(D)(l - D)} , 

i.e., it can be represented as a sum of independent contri- 
butions from different transverse modes with the trans- 
parencies D m . Similarly, the recurrence relation ( ^o|) and 
eq. (||) for the currents show that the same is true for the 
amplitudes A n . Therefore, the Fourier components (||) of 
the total current can be written as sums of independent 
contributions from individual transverse modes: 

h=Y J h{D m ), (12) 

m 

where the contribution of one (spin-degenerate) mode is: 
h{D) 



e 

itTi 



VDS ka - J detanh{^}(l - \a \ 2 )(a* 2k A* k - 



a- 2k A- k + ^(1 + a 2n a* {n+k) )(A n A* n+k - B n B* n 



(13) 



with the integral over e taken in large symmetric limits. 
The amplitudes B n in this equation are determined by 
the recurrence relation which follows directly from eq. 



n a2n+2Cl2n+l 

u— s a 



"2n+l 



a 2 n+l 
2 

2n+l 



D 



I — a. 

a2»Q2n^l p _ „l/2r 
: 5 O n -\ — —H OnO 



-)+l-4 n }B v 



R=l-D. (14) 
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FIG. 2. (a) DC IV characteristics and (b) differential con- 
ductance dl /dV of a short disordered SNS junction at several 
temperatures: T/A = 0, 1, 2, 3. (The T = 2A curve is omit- 
ted in (b).) The curves show the subharmonic gap singular- 
ities at eV = 2A/n, n — 1,2... associated with the multiple 
Andreev reflections. The differential conductance diverges as 
V -1 / 2 at small voltages - see text. 

Instead of using similar independent recurrence rela- 
tion for A n that follows from eq. (10), it is more conve- 
nient to determine these coefficients from an equivalent 
relation that can be obtained from a single-mode version 
ofeqs. (§) and (@) §|: 



A n +i — a,2n+ia<2nA n — 7? 1 / 2 (B 



(15) 
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The recurrence relations (|1J) and fllq ) can be solved 
explicitly in the case of perfect transmission, D = 1. In 
the general case of arbitrary D they provide an efficient 
algorithm for numerical evaluation of the current. There- 
fore, eqs. ([l^) - ( |l5| ) determine completely the time- 
dependent current in a short constriction with arbitrary 
distribution of transmission probabilities. In particular, 
we can use these equations to calculate the current in 
a short disordered SNS junction with large number of 
transverse modes N 3> 1 and diffusive electron transport 
in the N region. The distribution of transmission prob- 
abilities is then quasicontinuous, and is characterized by 
the density function p(D) (see, e.g., and references 
therein): 



dDp(D)... , p(D) 



1 



2e 2 £1(1 



L>)V2 ' 
(16) 



where G is the normal-state conductance of the N region. 
This distribution determines various transport properties 
of disordered mesoscopic conductors, for instance, the 
magnitude of the shot-noise suppression 



Figure 2 shows the results of the numerical calculations 
of the dc current-voltage (IV) characteristics and differ- 
ential conductance of the short SNS Josephson junction 
obtained by averaging the single-mode contributions to 
the current determined from eqs. (|l3|) - ( |l5| ) over the dis- 
tribution (|l|) . We see that the IV characteristics has all 
qualitative features of the high-transparency Josephson 
junctions: subgap current singularities at eV = 2A/n 
and excess current I ex at eV 3> 2A. It is instructive 
to compare quantitatively these features to those in the 
IV characteristics of a single-mode Josephson junctions 
plotted in Such a comparison shows that the mag- 
nitude of the excess current in the SNS junction, as well 
as the overall level of current in the sub-gap region cor- 
respond approximately to a single-mode junction with 
large transparency D ~ 0.8. At the same time, the sub- 
harmonic gap structure and the gap feature at eV = 2A 
are much more pronounced than in a single-mode junc- 
tion of this transparency. The amplitude of the oscilla- 
tions of the differential conductance corresponds roughly 
to the junction with D ~ 0.4 (although this comparison 
is not very accurate because of the different shapes of the 
curves). This "discrepancy" reflects the two-peak struc- 
ture of the transparency distribution ( |l6| ) of the diffusive 
conductor: the abundance of nearly ballistic modes leads 
to large excess and subgap currents, while the peak at 
low transparencies determines the SGS features. 

The main qualitative difference between the IV char- 
acteristics of the single-mode and diffusive SNS junction 
lies in the low-voltage (eV <§; A) behavior of the curves. 
As shown below, in contrast to a single-mode junction, 
the IV characteristics of the SNS junction have a square- 
root singularity at low voltages which directly reflects the 



shape of the high-transparency peak in the distribution 
dl6|). All this implies that the highly nonlinear character 
of the IV characteristics of the diffusive SNS junctions 
could serve as a direct experimental test of the trans- 
parency distribution (|l^) of the disordered mesoscopic 
conductor. 

Besides providing the basis for numerical evaluation of 
the current, the recurrence relations flli| ) and (|l5|) can be 
used to find the current analytically at small and large 
voltages. At large voltages (eV ^> A), the probability 
of MAR cycles decreases rapidly with the number of An- 
dreev reflections in them. We can then explicitly solve 
the recurrence relations limiting ouselves to the cycles 
with two Andreev reflections and find the current from 
eq. ([l3]) . For a single mode at T <C A we get the following 
high- voltage asymptote of the dc current pq] : 



D 2 



ln( 



i + Vr 



eD r AD, 

I(V) = — \eV+ (1- _ .... _,. 

K ' Trh 1 R 2VR(l+R) "l-VR 2eV 



)) 



(17) 



Averaging this result with the distribution ( |l6| ) we can 
obtain similar asymptote for an SNS junction: 

A 7T 2 A 2 

I(V)=G(V + -(--!)-—)}. (18) 
e 4 2eV 

The second term in this equation represents the excess 
current I ex and was first found in ]3^| by the quasiclassi- 
cal Green's function method. It can be checked that the 
asymptote (|l^) agrees well with the numerically calcu- 
lated zero-temperature IV characteristic shown in Fig. 
2a. The large-voltage approximation allows us also to 
find the asymptotes of the ac components of the current 
& 

The IV characteristics of the SNS junction can also be 
calculated analytically at small voltages, eV <C A, us- 
ing the understanding 23 3^] that the small finite volt- 
age V drives the Landau-Zener transitions between the 
Andreev-bound states of the modes with small reflection 
coefficients ]J« 1. For a single mode of this type, the 
nonequilibrium voltage-induced contribution to the cur- 
rent at T < A is §1: 



0, 

2exp{-?ri?A/e^} sin^/2, 



< If < it , 
7r < if < 2ir , 

(19) 



where f = fo + 2eVt/U is the Josephson phase difference 
between the superconductors. Averaging this equation 
with the distribution ( |l6[) and adding the equilibrium 
supercurrent from [p7| we obtain the dynamic current- 
phase relation of an SNS Josephson junction at eV <C A: 

/(^) = ^(co S (|)tanh- 1 [ S in(|)]+ 



4 



0. 



< (fi < TT . 



TT^/eV/Asm((f/2) , ir < Lp < 2ir . 



(20) 



The dc current / at small voltages is obtained by aver- 
aging eq. (|2^) over the phase (p: 



I = Gy/VA/e. 



(21) 



This square-root behavior of the current leads to the 
zero-bias singularity of the differential conductance of the 
SNS junction which can be seen in Fig. 2b. Physically, 
this large conductance is caused by overheating of elec- 
trons in the N region by the MAR process. Electrons 
with energies inside the energy gap traverse the constric- 
tion many times and as a result are accelerated to ener- 
gies much larger the eV. This means that the effective 
voltage drop across the constriction is much larger than 
V, leading to increased conductance. This mechanism of 
conductance enhancement is qualitatively similar to the 
so-called "stimulation of superconductivity" |58| (which 
is one of the plausible explanations of the zero-bias con- 
ductance singularities |39f| in long semiconductor Joseph- 
son junctions), although quantitatively the phenomena 
arc quite different. One of the most important differences 
between the short and long SNS junctions is that the sin- 
gularity ( pl| ) in short junctions should be suppressed by 
temperature simultaneously with the dc critical current, 
whereas in long junctions there is a temperature range 
where the zero-bias singularity is pronounced while the 
supercurrent is already negligible. 

The fact that the singularity (^lj) is caused by electron 
overheating implies that it should be regularized by any 
mechanism of inelastic scattering. Nevertheless, in con- 
strictions shorter than inelastic scattering length li n in 
the normal metal, the square-root conductance singular- 
ity should be experimentally observable at low tempera- 
tures, T <C A, when the inelastic scattering in supercon- 
ductors is also suppressed. 

All of the calculations presented in this Section were 
based on the assumption of ideal BCS electrodes of the 
constriction that are characterized by the Andreev reflec- 
tion amplitude ([!]). This assumption allows to describe 
dynamics of Andeev reflection with the Bogolyubov-de 
Gennes (BdG) equations. The most noticeable draw- 
back of this approach is the impossibility of incorporat- 
ing inelastic scattering in the electrodes which plays an 
important role in regularizing low-voltage singularities 
associated with MAR. In the next Section we show how 
the BdG approach can be generalized to superconductors 
with arbitrary microscopic structure which in general in- 
cludes inelastic processes. 



III. SUPERCONDUCTORS WITH GENERAL 
MICROSCOPIC STRUCTURE 

The method that allows to discuss the inelastic ef- 
fects in the superconducting electrodes of the constric- 
tion employs the quasiclassical non-equilibrium Green's 
functions. For the general introduction to this technique 
see, e.g., Jjoj . In this method, position dependent den- 
sity of states of the superconductors is described with the 
retarded and advanced Green's function Gr^a(x, e, e'). 
Deep inside the superconducting electrodes Gr^a reach 
their equilibrium position-independent form: 



9rM £ ) Jr,a( £ ) 
~fnA £ ) ~9rA £ ) 



S(e-s'), (22) 



where g and / are the normal and anomalous Green's 
function of the superconductor which satisfy the normal- 
ization condition: g 2 — f 2 = 1. Retarded and advanced 
functions are related simply as 

9A(e) = -g R (t) , /a(e) = -/fl(e) ■ 
Space dependence of Gr a is governed by the equation 

Mm 



isign(p x )v F 



dG 



R,A 



dx 



[Hr,a, Gr,a] 



(23) 



where vp is the Fermi velocity, sign(p x ) denotes two 
directions of propagation in the constriction, and H is 
the effective matrix Hamiltonian of the superconductor. 
Since equation ( p3| ) should be satisfied also in equilib- 
rium, H as a matrix commutes with G^ . Therefore, 
assuming electron-hole symmetry we can express it as 



Hr,a — ^R,aG^ a 



where f2(e) is complex quasiparticle energy. Imaginary 
part of fi comes from the imaginary part of electron self- 
energy and is positive: 



R,A 



> 0. 



(24) 



In order to calculate the current through the con- 
striction we need to find the Green's functions inside 
it. Equation ( p3| ) shows that the characteristic length 
scale of Green's function variation is the coherence length 
£ = hvp/A. In our model of short constriction with 
length d much smaller than £ this means that from the 
perspective of eq. ( p3| ) all points of the constriction cor- 
respond to x = 0. Equation ( p3| ) determines the Green's 
functions at x = through the condition that corrections 
G to the equilibrium functions (|2^) should decay inside 
the electrodes. To see how this condition translates into 
the matrix structure of G we perform a "rotation" in the 
electron-hole space which diagonalizes G^ (and hence 
H): 
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G R>A {e,e') -> U RtA (e)G R , A (e,e')U^(e') 
where the rotation matrix is 



(25) 



Ur,a = 



1 



1 " a R,A 



1 a R .A 
a R ,A 1 



Here 



Or(£) 



and = a^(e) . 



(26) 



We will see from the final results of this Section that a R 
defined by eq. ( poj ) has the meaning of the amplitude 
of Andreev reflection. In particular, for the ideal BCS 
superconductors the equilibrium Green's functions are: 
g R {e) = e/5, f R (e) = A/5, where 6 = [(e + zO) 2 - A 2 } 1 / 2 , 
and eq. ( p6| ) reduces to eq. of the previous Section. 

After the rotaion j25|) equation (^3[) for the retarded 
Green's function and p x > takes the form 



from the first row of the matrix equation ([28|) we obtain 
the recurrence relation for q\ : 

qi (e,e' + 2u) - qi (e,£')a R (s')a R (e' + u) = (5(e - e')+ 



+ g R (e))a R (e + u) - f R (e + u)a R {e')5(e - e' + 2u) . 

(29) 

Since the source terms in this relation are (5-functions we 
can look for a solution in the form qi — J2 n Ai(e)<5(£ — 
e' + 2un). The recurrence relation for the coefficients 
A n that follows from eq. (H|), complemented with the 
condition that A n should decay at n — > ±oo, can be 
solved directly. We obtain then the total Green's function 
in the constriction for p x > 0: 



Gr 



1 







-2a R {e) -1 



5(e 



2^5(8-8' 

n=l 



2un) 



(27) 



2n-l 



with tr's here and below denoting Pauli matrices. Com- 
bined with eq. (p4[), this equation means that the matrix 
structure of the solutions Gr\ }2 decaying in the first and 
second electrodes respectively, is: G\ri,2 °c a±, so that 



G m (e,e') = U R - 1 G m U R = «i(e,e') 



G fl2 (e,e') = U^G R2 U R = q 2 {e,e') 



a R (e') 1 
-a fl (£)aij(£') -Oii(e) 



-an(e) -a fl (e)ftfi(£^ 
1 Ofl(e') 



where q\^ 2 are some functions that will be found later. 

We discuss in detail the case of a single-mode ballistic 
constriction with D = 1. In this case the total Green's 
function should be continuous in the constriction: 



<-r R 2 + < ~ r R2 



(28) 



where G^j is the equilibrium Green's functions of the jth 
electrode. Expression (S) for these functions is valid 
only if the energies e, e' are measured relative to the 
Fermi energy of the corresponding electrode. If the zero 
of energy is chosen differently, the energies in eq. ( p^ ) 
should be shifted. For instance, if we chose the zero of 
energy to coincide with the Fermi level of the first elec- 
trode, then 



G 



(0) 

R2 



g R (8 + u)6(s - e') f R (8 + u)5(s-8' + 2u) 
-f R (e - u)S(e - e' - 2u) -g R {e - u)S(e ~ e') 



where u = eV , while G R \ is still given by eq. (|2^). Equa- 
tion ( p8| ) represents then the matrix equation that allows 
us to determine the functions qi, 2 . Indeed, eliminating q 2 



J\ a R (e+ju) 



3 = 1 



a R {e' 



-a R {e)a R (e') -a R {8) 



(30) 



In the same way we can find the Green's functions for 
backward propagation (p x < 0): 

Gr = ( I 2&K { e) ) S(e - e') + 2 % 2 ^) 



2n-l 

]J a R (e-ju) 

3 = 1 



clr(s) a R (s)a R (s') 
-1 -a R (e') 



(31) 



and also show explicitly that the advanced functions G a 
are related to G R as follows: 



G A (e, e', sign(p x )) = -G* r (e, e', -sign^)) 



(32) 



The current in the constriction depends not only on 
the density of states but also on the occupation of these 
states. The information about occupation probabilities 
is contained in the Keldysh Green's function Gk{e, s') 
which satisfy the following equation |4l],[l!| : 

QG 

isign(p x )v F -g— = HrPk + H K G A - G R H K - G K H A , 

(33) 

where Hra are the same matrices which determine the 
evolution of Gr,a in eq. (|23|), Hk = H R n — uHa, and 
n is the equilibrium quasiparticlc distribution, n(e, e') — 
'tanh(e/2T)(5(e- e'). 

Equation ( p3| ) implies that solution of eq. ( |33| ) can be 
written as 



Gk = G R n - uGa + G R , 



(34) 



G 



where Gh satisfies the homogeneous equation: 
dG 

isign(p x )v F — — = HrGh - GrHa- 
ox 



(35) 



Calculation of Gh follows closely the one for Gr,a that 
was described above. Diagonalizing Hr,a in eq. (|35| ) by 
transformation ( |25| ) and imposing the condition that Gh 
decays inside the electrodes, we see that for p x > and 
x = 0, Gri,2 oc (1 T a z) m the first and second electrode 
respectively. From this we find that 



Ghi(s,s') = Ur 1 GriUa 



1 Me') 
-Oij(£) -aij(e)aA(e') 



(36) 



h 2 (e,e') 



G H 2{e,e') = U r 1 G H 2Ua = 
-aR(s)a A {s') -a R (e) 



1 



oa(e') 



(37) 



Next we impose the continuity condition Gr\ = Gri 
in the constriction, and eliminating h 2 from this condi- 
tion obtain the recurrence relation for hi. This recur- 
rence relation can be solved similarly to the recurrence 
relation (|29|) and gives: 



h 1 (e,s') = 2j2 



2n+m 



5(e — e' + 2un) a R (e + ju) 



n=0 



3=1 



]]a A (e'+ju)n(-\e' 

3 = 1 



5 ( £ - e ' - 2un ) 



2n+m 



J|a fl (e + ju) } J a A {e' + ju)n ( \e + mu) 

3 = 1 3 = 1 



, (38) 



where > (e) = n(e + u) — n(e), and we used the conven- 
tion that II° =1 (...) = 1. Equations @, @ and @ 
determine the Keldysh Green's function in the constric- 
tion for p x > 0. Following the same steps we can also 
find Gr for p x < 0. 

The time-dependent current I(t) in the constriction 
can be expressed in terms of Gr as 

/(*) = g^/ & *' Tr h ( G ^ >0) ( £ ' e ')- 



-G^ <0) (e,e') 



exp 



i(e-n')t/h 



(39) 



For Gif found above (and given by the eqs. (p4|), (36), 
(|38|), ( ^0[ ) and (32)), eq. ( fj9| ) shows that the current is 
a combination of harmonics of the Josephson frequency 



ujj = 2eV/Ti, and the amplitude of the fcth harmonics is 
PI: 



e 



eVS h 



fcO 



oo m 



*tanh( — )(1- 



m+2k 



a.R{e) 



m— j — 1 j— rn+1 



(40) 



Transformations bringing expression for the current into 
this form used the relation o_r(— e) = —a* R (e). This rela- 
tion holds due to assumed electron-hole symmetry of the 
superconducting electrodes of the constriction. 
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FIG. 3. DC IV characteristics of a single-mode constric- 
tion between two superconductors with paramagnetic impu- 
rities for several values of (a) constriction transparency D, 
and (b) spin- flip scattering rate l/r s in the superconductors. 
G in normalization of the current is the normal-state conduc- 
tance of the constriction G = e 2 D/irh. The curves in (b) are 
shifted for clarity and correspond to (from bottom to top) 
h/r a A = 0.01, 0.1, 0.3, 1.0. For discussion see text. 



Comparison of eq. (gOJ) for the Fourier harmonics of the 
current in a single-mode constriction with the D — 1 ver- 
sion of eq. ( |l3| ) obtained from the BdG equations shows 
that these two equations coincide if we identify clr with 
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the amplitude of the Andreev reflection. To see this 
equivalence we note that at D = 1 the recurrence re- 
lations (|J) and @ give that B n = for all n, A n = 
for n < 0, and A n — 1 a i f or " > 0. With these 

wavefunction amplitudes eq. ( J13| ) indeed reduces to eq. 
(ff0|). This means that the only role of the complex in- 
ternal structure of the superconducting electrodes (in- 
cluding possible inelastic processes) in determining the 
ac Josephson current in a short ballistic constriction is 
to modify the amplitude of the Andreev reflection which 
is in general given by eq. (p6[). Starting from the Zait- 
sev's solution of the constriction problem in the matrix 
form p3] , one can show that this conclusion is valid 
for arbitrary D. The Fourier components of the time- 
dependent current in the constriction are determined by 
eq. (|l3|) and recurrence relations ([l4|) and (|l^ ) for general 
microscopic structure of the superconducting electrodes, 
if the amplitude of Andreev reflection cir is defined by 
eq. © @. 

To illustrate this approach, we calculated the IV char- 
acteristic of a short constriction between two supercon- 
ductors with paramagnetic impurities. In this case, the 
BCS singularity in the density of states at the gap edge 
is broadened by the spin-flip scattering and the magni- 
tude of the gap itself is suppressed. These two effects are 
reflected in the energy dependence of the Andreev re- 
flection amplitude au(e) which can be written in a form 
similar to eq. (Q): 



{e)=u{e)-{u(ef-lf/\ 



where u(e) is given by the equation 



e 
A 



h/r s A 



(41) 



(42) 



Here l/r s is the rate of the spin-flip scattering and A is 
the energy gap in absence of this scattering. Equations 
( [fl| ) and ( 42 ) allow for a simple numerical determination 
of an(e). DC IV characteristics for zero temperature and 
several scattering rates 1 /r s and transparencies D calcu- 
lated with the an(e) from eqs. (|il]) and ([l2j) are plotted in 
Fig. 3. Figure 3a shows how the IV characteristic evolves 
as a function of the constriction transparency D at weak 
spin-flip scattering. We see that even weak scattering re- 
sults in effective broadening of all current singularities. 
Due to suppression of the energy gap in the electrodes, 
the spin- flip scattering also gives rise to suppression of the 
zero-bias current jump at D = 1 and suppression of the 
excess current. Such a gap suppression becomes progres- 
sively more pronounced at larger scattering rates (Fig. 
3b), until the gap disappears completely at Ti/t s = A. 
In the gapless regime, the IV characteristic is practically 
linear. 



IV. NOISE 

The method developed in the previous Section allows 
us to calculate not only the average current in the con- 
striction, but also fluctuations of the current around the 
average. The spectral density of the current fluctuations 
Si(u>) in presence of the time-dependent average current 
is defined as 



1 

2^ 



dre lUT Ki(t,t + r), 



(43) 



K&t+r) = -(I(t)I(t+r)+I(t+T)I(t))-(I(t))(I(t+r)) , 

where the bar over Ki(t,t + r) denotes averaging over 
the time t. In this Section we discuss only the single- 
mode ballistic constriction. For such a constriction, the 
correlation function Kj can be expressed in terms of the 
quasiclassical Green's functions p5| : 

e 2 x 

A>(f,f + r) = --Y, Tr iG±(t,t + T)<j z Gi(t + r,t)a z + 



+G%(t,t + r)a,GX{t + r,t)a z ] , 

where ^2± is the sum over the two directions of propaga- 
tion in the constriction. The Green's functions G >:< are 
related to the functions Gr^a.k discussed in the previous 
section: 



1 



G >}< = t;[Gk±(G r -G a )} 



(44) 



The functions G(e, e') are the Fourier transforms of the 
G(t,t'). 

If we use in eq. (H) the Gr een's functions found in the 
previous section, we see that G >,<; can be written as 



(e, e') = J2 G «.'± (e + neV)5{e - e' + 2neV) 



(45) 



and the spectral density (M3) assumes the form: 



fdsTr[G> >± (e)a z G< n)± (s±fuj)a z 

" n, ±,±u> 



(46) 



Combinig this equation with the equations for the 
Green's function G >!< (eqs . (f44| ) and (E^) of this Section 
), (|30|) and (ft2Q of the previous 



and eqs. (|34|), Q3e 

Section) we arrive at the final result for Si (to) |2^ 
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[1 + 2ReJ2 II a ^ £ + 3 eV ) a A(z + jeV ± Hu)\ . (47) 
fe=li=i 

Function F here has the meaning of non-equilibrium dis- 
tribution of quasiparticles in the constriction: 

oo n 

F(s) = f(e) + E II I a «( £ ~ meV ) I' 

n—0 m— 

[f(s-(n + l)eV)- f(s-neV)], (48) 
where f(s) is the Fermi distribution function. 



fry=0.03A 



12. OO 



1=: 

5a 



8.00 




eV/A 

FIG. 4. Zero-frequency spectral density of fluctuations of 
the current in a single-mode ballistic constrictions at zero tem- 
perature as a function of the bias voltage. The peak at low 
voltages is due to multiple Andreev reflections. The inset 
shows a blowup of the peak at several temperatures: (from 
bottom to top) T/A = 0, 1, 2. 

The spectral density of current fluctuations defined by 
eqs. (|47]) and (|4^) depends strongly on the behavior of 
the Andreev reflection amplitude aa(s) as a function of 
energy e. Figure 4 shows the zero- frequency spectral den- 
sity S7(0) calculated numerically for the superconductors 
with small quasiparticle inelastic scattering rate 7, when 
an(s) is given by eq. ( p[i"| ) with 



u(e) 



e + ihj/2 



One can see from Fig. 4 that in this case the main feature 
of the current noise is the peak of 5/(0) at low voltages. 
The height of the peak is much larger than the classical 
shot noise result el /2ir which in our case would be on the 
order of e 2 A/27r 2 ?i. To understand the origin of this large 
noise we evaluate eqs. fl47| ) and ([48]) analytically in the 
limit V <C A/e. Expanding the amplitude a^(e) of An- 
dreev reflection in small relaxation rate 7 and replacing 
the sums with the integrals we obtain the spectral density 
of current fluctuations at low frequencies, lu <C A/h: 



SM = ^- r Y. [ A deF(e)(l-F(e±hu,)) 
n nV ± w J - A 



de' exp{- 



eWA 2 - v 2 



} 



cos(— [arccos(^) - arccos(^)]) . 



(49) 



where the quasiparticle distribution function reduces to 

dvfi'yiy) 



F(c) = /(e) 



(50) 



For ideal BCS superconductors 7^0 and eq. (|5C 
immediately gives that F(e) = /(A). Then the spectral 
density (Bfl) can be found explicitly: 



Si(lu) = 



1 + cos(t:Tiu) / eV) 



2^ 2 ?icosh 2 (A/2T)F (1 - {hiu/eV) 2 ) 2 



(51) 



This equation shows that the noise diverges at V — * 
as 1/V. At finite relaxation rate 7 this divergence sat- 
urates at V ~ h'y/e and gives the low- voltage peak in 
Fig. 4. The origin of this unusual voltage dependence of 
the noise is the process of multiple Andreev reflections. 
Each quasiparticle entering the constriction with the en- 
ergy equal to one of the gap edges generates an avalanche 
of Andreev reflections before it can escape out of the con- 
striction by climbing up or down in energy to the oppo- 
site edge of the energy gap. The number of generated 
Andreev reflections is 2A/ey, so that each quasiparticle 
causes a coherent transfer through the constriction of a 
charge quantum of magnitude 2A/V. For small voltages 
V <C A/e this charge is much larger than the charge of 
an individual Cooper pair. In this way the randomness 
of the quasiparticle scattering (quasiparticles get inside 
the energy gap with probability /(—A) from one elec- 
trode and with probability /(A) from the opposite elec- 
trode) is amplified. Therefore, the noise described by eq. 
( pl| ) can be interpreted as the shot noise of these large 
charge quanta, and we see that the process of multiple 
Andreev reflections dominates the current noise as well 
as the average current in short ballistic superconducting 
constrictions. 

This energy-domain picture with cycles of multiple An- 
dreev reflections, together with quantitative eqs. ( f49| ) and 
( pp| ) , has the "dual" time-domain formulation in terms of 
evolution of occupation probabilities of the two quasi- 
stationary Andreev-bound states localized in the con- 
striction that carry dc Josephson current Eq]. An ad- 
vantage of the time-domain formulation is that it can be 
generalized in a straightforward way to situations with 
the time-dependent bias voltage. It is also convenient 
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in establishing explicit relation between the ac, finite- 
voltage regime and dc regime in dynamics of supercon- 
ducting constrictions. In particular it shows how the 
finite- voltage current noise discussed in this Section goes 
over into large equilibrium supercurrent noise at vanish- 
ing bias voltage |2(||^J . 

V. CONCLUSIONS AND UNSOLVED PROBLEMS 

We have seen in the preseeding Sections that all elec- 
tron transport properties of short superconducting con- 
strictions are determined by the interplay of quasipar- 
ticle scattering inside the constriction and Andreev re- 
flection at the interfaces between the constriction and 
superconducting electrodes. Quantitatively, the scatter- 
ing process at finite bias voltages across the constriction 
is described by a set of recurrence relations for the am- 
plitudes of quasiparticle wave functions at energies e n 
shifted by integer number n of quanta of the Josephson 
oscillation, e n — s + 2eVn, n = 0,±1, .... The recurrence 
relations are valid for general scattering properties of the 
constriction that are characterized by the set of transmis- 
sion eigenvalues D m , and general microscopic structure 
of the superconducting electrodes characterized by the 
energy dependence of the Andreev reflection amplitude 
o(e). 

In this Chapter, we discussed the constrictions between 
identical s-wave superconductors. The recurrence rela- 
tions derived for this situation can be generalized to the 
case of two different superconductors (see the Chapter 
by P. Bagwell et al. in this volume) and d-wave super- 
conductors [ ^8p9| ]. An interesting possible direction for 
generalization of the recurrence relations is the case of 
transmission coefficient D(e) that is energy-dependent on 
the scale of superconducting energy gap A. Such a gener- 
alization would allow, for example, a systematic study of 
the ac Josephson effect through a resonant level J5(],|5l| or 
thick tunnel barrier. Another important unsolved prob- 
lem is the extension of the recurrence relations approach 
to the constrictions with the spin-dependent scattering, 
for example, Josephson junctions with ferromagnetic in- 
ter layer - see the contribution of G. Arnold et al. to this 
volume. 

The basic simplification which allows to describe quan- 
titatively such a variety of different situations is the ap- 
proximation of short constriction d <C £• A very impor- 
tant open problem is the lifting of this restriction. Gen- 
eral solution of this problem for junctions with d ~ £ ap- 
pears very difficult due to the necessity of self-consistent 
determination of the distribution of the pair potential 
A(a;) and electrostatic potential (f(x). As a first step to- 
wards such a solution, it would be interesting to apply 
the MAR approach to the opposite limit of long junctions 
d £. Long disordered SNS junctions have been attract- 
ing a lot of attention recently - see, e.g., and refer- 



ences therein, but the interest has been focused mainly on 
their properties close to equilibrium. MAR approach can 
be advantageous in studying strongly non-equilibrium 
situations. Another challenging unsolved problem is ap- 
plication of the results presented in this chapter to the 
description of Coulomb blockade in Josephson junctions 
with large transparency. The Coulomb blockade regime 
is characterized by quantum dynamics of the Josephson 
phase difference <p [|55|,^6| , in contrast to the classical dy- 
namics of tp assumed in this Chapter. 

In summary, one can expect further exciting new de- 
velopments and rapid progress in understanding of high- 
transparency Josephson junctions. 
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